clear;
hold off;
clearplot;
COTA = 0.000001;
MIN_ITER = 50;
rand("seed", 3001);

cant_numeros = 5000; %TODO: establecer criterio

%Genero la secuencia de numeros
i = 1;
do
	u(i) = rand();
	expu(i) = exp(u(i));
	
	if (i > 1)
		medias_u(i) = asimilar_media(medias_u(i - 1), u(i), i - 1) ;
		medias_expu(i) = asimilar_media(medias_expu(i - 1), expu(i),
			i - 1) ;
		covarianzas(i) = asimilar_covarianza(u, expu, medias_u(i),
			medias_expu(i));
		varianzas(i) = asimilar_varianza(mean(covarianzas), covarianzas,
	       		i);
	else
		medias_u(1) = u(1);
		medias_expu(1) = expu(1);
		covarianzas(1) = 0;
		varianzas(1) = 0;
	end
	i++;
until(i > MIN_ITER && abs(varianzas(i-1) - varianzas(i-2)) < COTA);

save ../datos/covarianzas.txt i covarianzas varianzas;
